An Effective Two-component Entanglement in Double-well Condensation 
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We propose a spin-half approximation method for two-component condensation in double wells to discuss 
the quantum entanglement of two components. This approximation is presented to be valid under stationary 
tunneling effect for odd particle number of each component. The evolution of the entanglement is found to be 
affected by the particle number both quantitatively and qualitatively. In detail, the maximal entanglement are 
shown to be hyperbolic like with respect to tunneling rate and time. To successively obtain large and long time 
sustained entanglement, the particle number should not be large. 
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I. INTRODUCTION 



The experimental achievement of atomic Bose-Einstein 
condensation(BEC) fT] has opened fascinating possibiHties 
for studying quantum properties of a macroscopic number of 
cold quantum atoms 1 2] and remains one of the most active re- 
search areas in recent years. Certainly, theoretical attentions 
are directed towards the underlying quantum correlation prop- 
erties of the condensed atoms |3]. Due to the macroscopic 
natural characteristics of EEC, it should be an ideal system 
for describing some quantum phenomena related to the co- 
herence. Correspondingly, although the Bose-Einstein con- 
densate is well represented by mean-field theory, it has many 
aspects that can be represented in a quantum picture contain- 
ing some proper description of correlations. The essential 
ultimate physics can be realized with the study of a simpler 
double-well EEC in the well-known quantum two-mode ap- 
proximation. 

As is well known, entanglement has come to be regarded 
as a physical resource which can be utilized to perform nu- 
merous tasks in quantum information processing. Also, apart 
from the fundamental physical interest in entanglement, the 
whole field of quantum computing and quantum information 
is based upon the ability to create and control entangled states 
|4J. In recent times, the study of the entanglement character- 
istics of various condensed matter system |5|] is focused on a 
pair of tunnel-coupled Eose-Einstein condensates (EEC's). In 
the simplest model, bosons are restricted to occupy one of two 
modes, each of which is in a EEC jQ]. 

A dynamical scheme of engineering many-particle entan- 
glement in EEC has been proposed by several authors, such 
as. Khan W. Mahmud et at They introduced a quantum 
phase-space dynamics to generate tunable entangled number 
states using the Husimi projection method |8]. However, the 
notion of entanglement in macroscopic ensembles allows to 
investigate the boundary between quantum physics and clas- 
sical physics and, possibly, could also give some insight into 
the measurement process |9, 10, 11]. In this paper we will 
investigate in detail a scheme that measures entanglement of 
N two-component particles trapped in a double- well under a 




FIG. 1 : A sketch of the double-well system with two-component condensate 
trapped in. Initially, the atoms of the component A are trapped in the left side 
and the atoms of the component B are trapped in the right side. 



kind of reasonable approximate assumption. This approxima- 
tion is presented to be valid under stationary tunneling effect 
for odd particle number of each component. 



II. THE MODEL 

Our system consists of a double-well and two component 
condensate trapped in it as depicted in FIG. 1. Initially, the 
atoms of component A and E are separately located in left 
and right potential wells, respectively, fT?]. The many -body 
Hamiltonian for a system of weakly interacting bosons in 
an external potential V{r), in second quantization, is given by 
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where i^(r) and i^'(r) are the boson annihilation and cre- 
ation field operators, m is the particle mass and g - '^""^f , 
with a.5 the s-wave scattering length. In studies of double- 
well EEC or two-component condensates, the low-energy 
many -body Hamiltonian in Eq. (I) can be simplified in 
the well-known two-mode approximation fisl [3, which 
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has been used widely in studying the double-well conden- 
sate. Under this approximation, the system is modeled by the 
Hamiltonian(/i = 1) 112], 

H = ^(fl^flfi -I- fl^«L) + ^(.K^R + b^Ji) + Kialaiblbi 

+ Y[(btbLf + (KbR)\ (2) 

where the subscripts L and R denote respectively the local- 
ized modes in the left and right potential wells. Because of 
two modes available for each component, there are four oper- 
ators fl| and (j - L, R) that denote the creation operators 
of the components A and B in two wells in the model. The 
parameters O,-, /c,- with / = A,B and k describe the tunneling 
rate, self-interaction strength of the component A (B) and the 
interspecies interaction strength. Especially, Q, depends lin- 
early on particle number A^,, so that, Q, = £lr + gTi(Ni - 1), 
where eiu, T\ are the parameters after second quantization and 
A^, is particle number of component ; J3l • For simplicity, we 
set = = Q in following sections. 



III. CLASSIC ANALYSIS 

Due to the macroscopic nature of its wave function, BEC 
should be an ideal system for the generation of macroscopic 
quantum superposition states (Schrodinger cat states) 131. 
Such a system could be, on one hand, analyzed in a classi- 
cal method from a large particle number point of view and on 
the other hand studied in a quantum method because of the 
quantum superposition character. 

The classical Hamiltonian that describes the mean-field dy- 
namics of BEC in a double well has been studied in sev- 
eral papers In a mean-field assumption for the 
two-mode double well BEC in case of large enough A^, the 
operators cij can be replaced by a c-number -^nje'"', where 
j = L,R. With this assumption and defining population dif- 
ference (POD) n = ^^^^2^, = - 0R, we analogously 
introduce c-numbers both for component A and B, that is. 



nAL + nAR = Na, 

nBL + nBR = Nb, 

nAL - nAR 
nA = r , 



riB - 



nBL - nBR 



the classical Hamiltonian is then given by 



^Nl - {Iha)^ cos 6 a + ^^Nl - ilnB^ cos Ob 



K Ki -, Kb ^ 

+-NaNb + -fNl + -fNl 

+KAn\ + KBn\ + iKHAns 




FIG. 2: Shown is the population difference for component A with parameter 
values K = 20, ka = 20, Kg = 20. 





FIG. 3: Shown is the population difference for component B for parameter 
values *■ = 20, ka = 20, kb = 20. 



and the equations of motion are 



nA 
iiB 



= -^Nl-{2nAf sin^A, 
= j^Nl-ilriB)^ sm0B, 



IDiiA COS 6a 
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+ iKAnA + 2m B, 



Nl - (2nA)^ 
2QnBCOs 6b 



4 



+ 2KBnB + 2KnA 



(4) 
(5) 
(6) 

(7) 



A^^ - (2nB)2 



(3) 



We plot the amount of the POD for component A and B 
{nA,nB) as a function of the tunneling rate Q. and time t in 
FIG. 2 and FIG. 3. 

FIG. 2 shows the POD for component A with parameter 
values K = 20, ka = 20, kb = 20. The POD between two wells 
of each component presents a behavior of local maximum and 
minimum. The period is determined by the particle number 
of the component. The larger the particle number, the smaller 
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the period. It should be stressed that the peak of the POD of 
component A between two wells is decreased gradually. In 
FIG. 3, there have the same conditions for component B but 
for the negative initial value of POD. 

FIG. 2 suggests that the stationary state of population would 
be in particle number equilibrium. Similar result has been pre- 
sented in Ref. [12], where the expectation value of the defined 
inversion operator {a^ai - o^or} was found to tend to zero. 
Our solution is based on a continuous analysis of Equ. [4- 
7]. While, in physical system scenario, the population number 
must be discrete since there are 1 orthogonal Fork states of 
each component. That is, for even component particle number 
A^, the stationary population in each well would be j. While 
for odd A^, the stationary POD would be ± j . Whatever be the 
initial population number of components in each well, a min- 
imal POD-0 or 1 can be obtained due to the tunneling effect 
and repulsive interaction. 

In the following section, we shall see how to generate and 
measure a mode entangled state under such an approximate 
situation. 



IV. GENERATING AND MEASURING OF 
ENTANGLEMENT UNDER SPIN-HALF APPROXIMATION 
ASSUMPTION 

The quantum phase-space model presented in Ref. [7] 
points to a way that an entangled state can be generated with 
a single-component BEC in double well. The authors showed 
the generation of tunable entangled states in phase space us- 
ing Husimi distribution function. We now introduce another 
method to prepare entangled states. Schwinger has developed 
an entire angular- momentum algebra in terms of two sets (up 
and down) of creation and annihilation operators for uncorre- 
cted harmonic oscillator constructing the angular-momentum 
operators as fnll 

J+ — J.x + iJv — a^a^, 
J- - Jx - iJy - a^a+, 
1 , 

Jz = 2^a_^_a+ - a_a-). (8) 

This construction then satisfies the standard angular- 
momentum commutation relation [Ji,J„,] - iei„,„J„. Now, we 
regard the double-well as a two-state system and similarly in- 
troduce an analogous-angular momentum algebra to express 
the Hamiltonian in Eq. [2]. For the two-component conden- 
sate trapped in a double well,we defining 

5r = ALA^, 



S\ = -(A^Az, - A^A«), 



(9) 



where /=A, 
tively). 5 + 



B, A = a,b (for components A and B respec- 
5^ satisfy the angular-momentum commutation 

relations: 1^5^, = 25^. With these descriptions, we, in 
fact, consider the left and right wells as two states (|L) and 
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FIG. 4: Shown is the distribution of pseudo-spin 5; with the respect to POD. 
Solid line for odd particle number, dotted line for even particle number. 



\R)) of the particles. As for A^-particle system, the collective 
spin is S, = ~f ■ ■ ■ f - "^o, there would be A^,- "excited states" 
for component / with particle number A^, as illustrated in FIG. 
4. Now, we assume the initial particle number of each compo- 
nent is odd. The stationary POD, from the analysis in the last 
section, must be 1 . So, the collective spin of each component 
should be This assumption suggests us to eliminate adi- 
abatically "high energy" states of each component, thus each 
component can be treated as a pseudo-particle with expecta- 
tion value + i of spin. 

Then, the raising and lowering operators of this pseudo- 
particle 5^ (5 ) can be expressed in a 2 x 2 dimensional 
Hilbert space by 



\l}a{RI 

\L)b{R\, 



\R}a{li 

\R)b{L\- 



(10) 



After the operators at (at) replaced by the S^{Sj), where 
j - A, B, i - L, R, the global Hamiltonian in Equ. [2] can be 
spanned on the basis of \L)a ® \L)b, \L}a ® \R)b, \R)a ® \L)b, 
\R)a ® \R}b- We obtain the efifective Hamiltonian of our system 
as 



H ■■ 



( Ki K K ^ 

K K2 K 

K Q Ki K 

K K Ki J 



(11) 



where K ^ 'j, Ki ^ k + !^ + K2 ^ 



2 . The eigenvalues 
and eigenvectors of the Hamiltonian can be found explicitly. 



El ^K2, <pi ^—^i\L)A\R)B- 
V2 



\R)a\L)b); 



£2,3 



02,3 - ^2,; 



K2 + KI + & 

2 ' 
\L)a\L)b^ 



-K2 + Ki-e 

4K 



i\L)A\R)B + \R)a\L)b) + \R)a\R)b 



£4 = ^^1, <p4 = -^i\L)A\L)B - \R)a\R)b). 
V2 



(12) 
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FIG. 5 : Shown is the degree of entanglement as a function of the tunneling 
rate Q. and time t under the two-mode approximate Hamiltonian for parameter 
values K = 20, /f^ = 20, Kg = 20. 




Q. 



FIG. 6: Shown is the maximal entanglement for different k as a function of 
the tunneling rate SI and time f under the two-mode approximate Hamiltonian 
for parameter values k = 20, ka = 20, kb = 20. 



where ^2, ^3 are normalized factors and 
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Kl - 2K1K2 + K\ + \6K^. We assume that, initially, 

the POD of component A is |, and the POD of component B 
is -i, so that, the initial state of system is 



l^^(O)) = \L)a ® \R)b. 



(13) 



Wootters Concurrence has been widely used in measuring 
the entanglement of bipartite two-state system, which is de- 
fined as L18.1 . 



C(p) - max {Q, A\ - A2 - At, - /I4}, 



(14) 




where /I,- are the square roots of he eigenvalues of the non- 
Hermitian matrix pp with p = (cr,, ® cry)p*(cry ® cr,,) in de- 



FIG. 7: Shown is the degree of entanglement as a function of time t under 
the two-mode approximate Hamiltonian for parameter values n = I, k = 20, 
= 20, Ks = 20. 



creasing order. Wootters Concurrence gives an explicit ex- 
pression for the entanglement of formation, which quantifies 
the resources needed to create entangled state. 

Then the behavior of entanglement in this system described 
by spin-half assumption Hamiltonian is illustrated in FIG. 5 
and FIG. 6. In this situation when there is no decoherence, 
entangled states can be generated because of the overlap of 
the wave packet of each component in two wells. We note that 
the entanglement presents a periodic distribution for any value 
of t falling in the region, that is, the behavior of the amount 
of entanglement between the two modes is non-mono tonic, 
the maximum entanglement can be achieved if we control the 
tunneling time at about t - -1- y). The corresponding 
maximal entanglement for different k described by FIG. 6 is 
shown to be hyperbolic like with respect to tunneling rate and 
time. The half width of the peak of entanglement t can be 
defined as a parameter to represent the quantity and the quality 
of the maximal entanglement. As an example, we show the 
entanglement when Q = 1 in FIG. 6. At the peak position 
of the entanglement, t « + y), the half width is t a 

31.369. While for Q = 1.5, the half width is t * 13.899, 
which suggests that the larger the Q, the smaller the half width 
of the peak. Then, we can easily obtain large and long time 
sustained entanglement successively for 1 < Q < 1.5, which 
presents a broad banding distribution. 



V. CONCLUDING REMARKS 

The creation of many particle entangled states in macro- 
scopic systems is one of the major goals in the studies on 
fundamental aspects of quantum theory. In this paper, the 
studies we have performed show a spin-half approximation 
method for two-component condensation in double wells to 
discuss the quantum entanglement of two components, which 
is presented to be valid under stationary tunneling effect for 
odd particle number of each component. The evolution of the 
entanglement is found to be affected by the particle number 
both quantitatively and qualitatively. In detail, the maximal 
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entanglement are shown to be hyperbolic like with respect to 
tunneling rate and time. To successively obtain large and long 
time sustained entanglement, the particle number should not 
be large. Thus we can transform an exact many-body prob- 
lem which is difficult even for single-component condensate 
to a bipartite two-state problem, similarly, we could measure 
the entanglement simply and conveniently. There also refer 
to the comparison of classical and quantum dynamics, great 
progress has been made in this subject i^, which is inter- 
esting and sensitive recently and is good for further studies. 
Our entanglement states are specified with the particle spatial 
location in left or right well rather than internal energy lev- 



els, this may be operational in physical applications such as 
quantum entangled particles distribution and quantum mea- 
surement. Moreover, the method developed here may find ap- 
plications in the studies of entanglement of other BEC. 
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